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ABSTRACT 


Studies of the GPS satellite orbit determination problem which have been 
published previously have assumed that the tracking data would he provided by 
pseudo-range measurements i:rom about 4 ground stations. Substantial allow- 
ances were Included in the error budget for drifts in the satellite and ground 
station clocks, and for uncertainty in the water vapor contribution to the 
tropospheric propagation correction. We have started an investigation re- 
cently to see how much improvement in the orbit accuracy can be expected if 
hlgh-accuracy tracking data from a substantially larger number of ground sta- 
tions is available. For observations from 20 ground stations, the results 
Indicate that 20 cm or better accuracy can be achieved for the horizontal 
coordinates of the GPS satellites. With this accuracy, the contribution to 
the error budget for determining 1000 km baselines by GPS geodetic receivers 
would be only about I cm. 


INTRODUCTION 


A preliminary study of the accuracy with which geodetic baselines can be 
determined using signals from the Global Positioning System (GPS) satellites 
was very encouraging (Bossier et al . 1980), and we have since carried out a 
much more detailed study (Larden and Bender, 1982a). It was assumed that the 
GPS receivers would measure the phase of the reconstructed carrier signals 
which are generated in some receivers by using the modulation code (Ward, 
1982) and in other instruments without requiring knowledge of the code 
(MacDoren et al . , 1982; Counselman and Steinbrecher, 1982). Baselines up to 
about too km long were considered, and uncertalnites in the horizontal 
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coordlaates £or the GPS satcillccs were taken to be 2 meters. The measurement 
accuracy was assumed to bo 1 cm for the range difference from two ground sta- 
tions to a given satellite, with most of the error due to the expected uncer- 
tainty in the tropospheric correction even if water vapor radiometers are used 
at both ground stations to determine the integrated water vapor content along 
the line of sight. A “modified worst case" type of analysis was employed be- 
cause of the expected substantial correlation of the tropospheric correction 
errors over the measurement times. 

With the above assumptions and analysis method, the expected baseline 
uncertainties were found to be less than 1 cm for horizontal components and 2 
cm for the vertical for short baselines in most parts of the world, even with 
only the initial 6 GPS satellites in operation. However, for baselines sub- 
stantially longer than 100 km, the largest part of the uncertainty would come 
from the assumed GPS satellite position uncertainty of 2 meters for the hori- 
zontal coordinates. This was one reason why we decided to investigate the 
extent to which determination of the GPS satellite orbits could be Improved if 
accurate reconstructed carrier phase measurements were made at a number of 
worldwide tracking sites whose coordinates wer" known from other types of 
measurements (harden and Bender, 1981). The other ’■eason was the probable 
major Importance of receiving GPS signals aboard altimeter satellites such as 
the TOPEX (Topographic Experiment) satellite and the ERS (Earth Resources 
Satelllta) in order to determine their orbits accurately, provided that the 
GPS satellite locations could be determined accurately (Bender and harden, 
1982). The initial results of our study of the GPS orbit determination pro- 
blem are given in the remainder of this paper. 

GPS SATEhhITE ORBIT DETERMINATION 


The objective of our study is to establish the accuracy achievable for 
determining the GPS satellite orbits when all 18 of the planned satellites are 
in operation. IF the GPS signals are as useful os expected for geodetic 
applications and for determining the orbits of other satellites, then the 
tracking of the GPS satellites from as many as 20 fixed ground stations seems 
likely. Such a future network might be a combination and extension of the 
present TRANET 2 and MEDOC networks, with each station equipped with a high- 
accuracy GPS receiver and a water vapor radiometer. 

Since the goals of the NASA Crustal Dynamics Project and of related 
projects in other countries include determining intercontinental baselines 
with 3 cm accuracy, we assume 3 cm uncertainty in each coordinate for the GPS 
tracking stations. It appears quite feasible to achieve this accuracy with 
both LAGEOS ranging and VLBI measurements using highly mobile stations. The 
additional cost of such measurements appears t6 be reasonable, in view of the 
expected worldwide benefits, provided that the main part of the Crustal 
Dynamics Project or some ::imllar program is being carried out. Repeat mea- 
surements would be needed only every few years for most stations, provided 
that the relative plate motions are being monitored and that the expected low 
rates of intraplate motion are found in most areas. 

The main problem in determining the expected accuracy of the GPS orbits 
is chat Che uost serious limitation is likely to come from time-dependent 
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systematic errors whose time variations are not known and cannot be modeled. 

In order to stake some appropriate allowance for such errors, we assume that 
measurements are made at each ground station only once every half hour. For 
orbital arcs of 2 and 4 hours, this means that either 3 or 9 sets of measure- 
ments are Included In the simulations. The measurements acutally would be 
made much more frequently, and the part of the error which is truly uncor- 
relatcd from measurement to measurement would nearly average out. Thus the 
assumed random measurement errors really are intended to represent the effects 
of the time-varying systematic errors. 

In our simulations we assume that reconstructed carrier phase measure- 
ments are made at nearly the same time from each ground station to all GPS 
satellites which are above 20“ elevation angle. The time spread for the dif- 
ferent measurements with respect to the nominal measurement tlsie is assumed to 
be short enough so that only a single clock correction needs to be considered 
for each satellite ami each ground station. We take the difference In range 
from a particular satellite to a given pair of ground stations as the obser- 
vable, so that corrections to the satelllre clocks do not appear explicitly In 
the analysis. Corrections to the clock differences between 19 of the ground 
station clocks and one reference clock are solved for at each nominal mea- 
surement time. The only other parameters solved for are the six orbit ele- 
ments for each satellite. 

The measurement noise was taken to be 1.4 cm for the difference In range 
to a given satellite between a reference station and each of the other sta- 
tions for which the elevation angle was greater than 20“ . The station with 
the highest elevation angle for the satellite usually was chosen as the 
reference station. Other error sources allowed for were a 5Z uncertainty in 
thCgSolar radiation pressure force acting radially away from the sun, a 5 x 
10 fractional uncertainty In Wl, where M Is the earth's mass, and the GEM-10 
scaled standard deviations for each of the gravitational field harmonic coef- 
ficients through degree and order 4. An additional thrust uncertainty equal 
In magnitude to the 52 uncertainty in the solar radiation pressure force 
usually was included, with Its direction making equal angles with the along- 
track, cross-track, and radial directions. 

The 20 ground station locations chosen for use In the simulations are 
listed in Table I. The first 3 are currently GPS monitor station sites, and 
the next 6 are TRANET sites. The following 3, Bermuda, Reunion, and Unalaska 
are close enough to the TRANET sites In Herdon, Virginia, In Mahe, Seychelles 
Islands, and In Anchorage, Alaska so that these sites probably would give 
comparable accuracy. The following 2 in Djibouti and Pretoria are MEOOC 
sites, and the remaining six are not presently connected with Transit or GPS 
satellite tracking. Auckland and Geraldton are Included Instead of the TRANET 
site In Adelaide, Australia, in order to give better geometrical coverage. 
Dakar, New Delhi, Quito, and Rio Gallegos represent general areas where 
coverage vAould be desirable. 
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Table 1. 

Assumed GPS cracking station locations 


I < Guam 
2. Hawaii 

3« Vandenberg AFB, Calif* 

4 . Grasse , France 

5. Kushiro, Japan 

6. McMurdo Sound, Antarctic 

7. Pretoria, South Africa 

8. Sao Paulo, Brazil 

9. Thule, Greenland 

10 . Bermuda Islands 


11. Reunion Island 

12. Unalaska, Aleutian Islands 

13 Djibouti, Republic of Djibouti 

14. Papeete, Tahiti 

15. Auckland, New Zealand 

16. Dakar, Senegal 

17. Geraldton, Australia 

18. New Delhi, India 

19. Quito, Equador 

20. Rio Gallegos, Argentina 


RESULTS AND CONCLUSIONS 

The results obtained from the main simulations are given in Table 2. The 
uncertainties in the GPS satellite horizontal coordinates are the values cor- 
responding to the middles of the arcs. It can be seen Chat the dominant 
errors are those due to the unmodeled parameters. Of these, the major part 
comes from the assumed 3 cm vmcertalnty in each component of the >;round 
station locations. The uncertainties in GM and in the otner low degree 
gravity field coefficients give a total of only roughly 3 cm uncertainty in 
each of the horizontal coordinates in the 4 hr case and abovit half as much in 
the 2 hr case. The uncertainties due to solar radiation pressure and to the 
additional assumed thrusts also are small. The contributions from the random 
errors were scaled up by a factor (k/4)^' , where k is the number of 
measurement times. This corresponds to normalizing the results to the rather 
arbitrarily chosen case of 4 observation times, so that correlated errors 
don't tend to average out as much. 

From the results given, it appears that uncertainties of less than 20 cm 
for the horizontal coordinates of the GPS satellites probably can be achieved 
using the 4 hr arcs, provided that a sufficient number of ground stations 
equipped with the right kinds of receivers and water vapor milometers provide 
regular tracking data. The uncertainties obtained for the 2 hr arcs were 
about 35% larger. The model assumed for radiation pressure force uncertain- 
ties was not realistic enough to give good estimates of the effect of this 
error source, since the actual error in the satellite acceleration will not be 
confined to the two directions included in the model and will vary appreciably 
over the 4 hr arc. However, the magnitudes of the corresponding orbit errors 
were small enough so that the radiation pressure effects do not appear to be 
serious limitations on the accuracy for arcs of this length. 

In view of the results obtained, it appears that GPS receivers can be 
used to obtain accurate positions over baselines up to at least 1000 km in 
length. For this distance, the contributions of our horizontal coordinate 
uncertainties of less than 20 cm for the GPS satellites to the baseline 
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Table 2. 

Uncertilntf.es In GPL satellite horizontal coordinates 


Contributions to 



Position Uncertainties (cm) 


total uncertainty 


2 hour 

** 

case 



4 hour 

case 

■ 

Cross 

track 

Along 

track 

Cross 

track 

Along track 

Modified random 
noise,* RMS 

10.5 


11.1 


6.4 


6.6 

Unraodeled effects, 
RMS 

18.0 


19.1 


14.9 


14.3 

RSS total 

20.8 


22.1 


16.2 


15.7 

Worst Individual 
uncertainty for 
any satellite 

25.3 


26.4 


18.5 


18.9 


^ 1/2 
Random noise contributions were multiplied by a factor (k/4) , where k is 

^j^he number of measurement times* 

The thrust uncertainties were not included In this case. 


uncertainty can be shown to be only about 1 cm In each coordinate (harden and 
Bender, 1982b). Thus It Is likely to be useful for a number of countries to 
participate In accurate tracking of the GPS satellites, so that the distances 
over which GPS receivers can be used for hlgh-accuracy geophysical and geo-* 
detlc measurements can be considerably extended. The use of Improved GPS 
satellite positions for determining the orbits of the TOPEX and ERS satellites 
and other altimeter satellites also is likely to be of major importance 
(Bender and harden, 1982; harden and Bender, 1982b). 
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